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Abstract: We propose a renormalizable model of a fermionic dark matter by introducing 
a gauge singlet Dirac fermion and a real singlet scalar. The bridges between the singlet 
sector and the standard model sector are only the singlet scalar interaction terms with the 
standard model Higgs field. The singlet fermion couples to the standard model particles 
through the mixing between the standard model Higgs and singlet scalar and is naturally 
a weakly interacting massive particle (WIMP). The measured relic abundance can be ex- 
plained by the singlet fermionic dark matter as the WIMP within this model. Collider 
implication of the singlet fermionic dark matter is also discussed. Predicted is the elastic 
scattering cross section of the singlet fermion into target nuclei for a direct detection of the 
dark matter. Search of the direct detection of the dark matter provides severe constraints 
on the parameters of our model. 
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1. Introduction 

The missing mass of some non-visible form of matter in the galaxy cluster was first inves- 
tigated by Zwicky in 1933 JTJ. Since then, there have been a lot of efforts to probe the 
dark matter (DM) which provides the unseen mass in the cluster. Evidences have been 
found including the galactic rotation curve Q and the observation of the Bullet cluster || . 
The precise measurement of the relic abundance of the cold dark matter (CDM) has been 
obtained from the Wilkinson microwave anisotropy probe (WMAP) data on the cosmic 
microwave background radiation as Q 

0.085 < n CD Aih 2 < 0.119, (2a level) (1.1) 

where is the normalized relic density and the scaled Hubble constant h ~ 0.7 in the units 
of 100 km/sec/Mpc. 

Since there is no proper CDM candidate in the standard model (SM) contents, the 
extended models of the SM are required to provide a candidate of DM. Various candidates 
of the CDM have been proposed. Weakly interacting massive particles (WIMP) are favored 
to explain the observed value of the relic abundance in view of new physics beyond the SM. 
WIMPs include the lightest supersymmetric particle (LSP) in the supersymmetric models 
with R parity (B|, [|, the lightest Kaluza-Klein particle in the extra dimensional models 
with conserved KK parity 0, and the lightest T-odd particle in the T-parity conserved 
little Higgs model ||. Addition of a real singlet scalar field to the SM with Z2-parity 
has been considered as one of the simplest extensions of the SM with the nonbaryonic 
CDM H 1C, 11]. A model which introduces singlet Majorana neutrinos and a singlet 



scalar has been considered in Ref. [12]. Another scalar extension with multiple Higgs 



doublets is presented in Ref. [13|. General classification of the extra gauge multiplets as 



a minimal dark matter candidate has been performed in Ref. [14]. On the other hand 



incorporating the gauge coupling unification with the dark matter issue, a fermionic DM 



with the quantum numbers of SUSY higgsinos and a singlet is suggested [15, 16]. Heavy 
particle DM communicating with supersymmetric SM via pure Higgs sector interaction, 
which is motivated by Higgs portal and Hidden valley models is also considered] 17 1. A 
model with a gauge singlet Dirac fermion is proposed as a minimal model of fermionic dark 
matter p!^]. In this model, the singlet fermion interacts with the SM sector only through 
nonrenormalizable interactions among which the leading interaction term is given by the 
dimension five term (1/K)H^ H-tptp, where H is the SM Higgs doublet and tp is the dark 
matter fermion, suppressed by a new physics scale A. 

We propose a renormalizable extension of the SM with a hidden sector which consists 
of SM gauge singlets (a singlet scalar and a singlet Dirac fermion). The singlet scalar 
interacts with the SM sector through the triple and quartic scalar interactions. There are 
no renormalizable interaction terms between the singlet fermion and the SM particles but 
the singlet fermion interacts with the SM matters only via the singlet scalar. Therefore it 
is natural that the singlet fermion is a WIMP and a candidate of the CDM. Our model is a 
minimal model of renormalizable extension of the SM including the fermionic dark matter. 
The model is described in section We show that the singlet fermion can be a CDM 
candidate, which explains the measured relic density by the WMAP with the experimental 
constraint on the Higgs bosons at LEP2 in section || The direct detection of the fermionic 
CDM is investigated in Sec.||. Finally we conclude in Sec.||. 



2. The model 

We introduce a hidden sector consisting of a real scalar field S and a Dirac fermion field 
-0 which are SM gauge singlets. The singlet scalar S couples to the SM particles only 
through triple and quartic terms with the SM Higgs boson such as SH^H and S 2 H^ H. 
New fermion number of the singlet fermion is required to be conserved in order to avoid 
the mixing between the singlet fermion and the SM fermions. The global U(l) charge 
of the singlet Dirac fermion takes the role of the new fermion number. As a result, no 
renormalizable interaction terms between the singlet fermion ip and the SM particles are 
allowed. Thus the interaction of tp with the SM particles just comes via the singlet scalar. 
We write the Lagrangian as 

£ = C-SM + £hid + £%nU (2-1) 

where the hidden sector Lagrangian is given by 

£hid = Cs + Cj,- gs4>ipS, (2.2) 

with 

c s = \ M {&>S) - ^S 2 - ^ - 

= i>(i<j) -m^ )i/j. (2.3) 



The interaction Lagrangian between the hidden sector and the SM fields is given by 

C mt = -\ X H^HS - X 2 H^HS 2 . (2.4) 



The scalar potential given in Eq. fl2.3| ) and ( |2,4[) together with the SM Higgs potential 
—[i 2 H^H + Xq(H^ H) 2 derives the vacuum expectation values (VEVs) 



(H) 




(2.5) 



for the SM Higgs doublet to give rise to the electroweak symmetry breaking, and (5) = 
xq for the singlet scalar sector. The extremum conditions 9V/dH\ <H0>=VQ ,^2 = and 
dV / dS\ < s>=x = lead to the relations |H| 



H 2 = \ v 2 + (Ai + X 2 x Q )xo, 
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(2.6) 



The neutral scalar states h and s defined by H° = (vq + h) / \[2 and S = xq + s are mixed 
to yield the mass matrix given by 



2 = zX. 
Mft ~ dh 2 

2= &v 

Ms - ds 2 
d 2 V 



h=s=0 



As = 



dhds 



A 3 A 4 2 Ai^o 

XQ H Xn . 

2 3 2x 

= (Ai + 2X 2 xo)v . 



(2.7) 



/i=s=0 



The mass eigenstates h± and /12 are obtained by 

h\ = sin 9 s + cos /i, 
/12 = cos # s — sin /i, 

where the mixing angle 9 is defined by 

tan 9 



1 + ^TTy 2 

with y = 2/j| s /(/j| — fi 2 ). The Higgs boson masses mi and m,2 are given by 



(2.8) 



(2.9) 
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where the upper (lower) sign corresponds to mi (7712). According to the definition of tan#, 
we get I cos 9\ > implying that hi is SM Higgs-like while hi is the singlet-like scalars. As 
a result, there exist two neutral Higgs bosons in our model and the collider phenomenology 
of the Higgs sector might be affected. We will discuss it in later section. 

The singlet fermion tp has the mass mj, = m^ + gs%o as an independent parameter 
of the model since m^ is just a free parameter. The Yukawa coupling g$ measures the 
interaction of ip with other particles. Generically the interactions between ijj and the SM 
particles are suppressed by the mass of singlet scalar and/or the Higgs mixing. Therefore 
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if? is naturally weakly interacting and can play the role of a cold dark matter as an WIMP. 
If we fix masses of two Higgs bosons, the singlet fermion annihilation processes into the SM 
particles depend upon the fermion mass Yukawa coupling gs, and the Higgs mixing 
angle 0. If the final state includes Higgs bosons, h\ or /12, several Higgs self-couplings are 
involved depending on various couplings in the scalar potential. 

3. Implications in cosmology and collider physics 



In the early universe, the dark matter is assumed to be in thermal equilibrium by the active 
annihilation and production process with the SM sector. When the universe cools down 
and the temperature T drops below the DM mass, the DM number density is suppressed 
exponentially so that the annihilation rate of the dark matter becomes smaller than even 
the Hubble parameter. Then the interactions of the DM freeze out to make the DM 
particles fall out of equilibrium and the DM number density in a comoving volume remains 
constant. Therefore the current relic abundance of the CDM depends on the annihilation 
cross section of ip into the SM particles or the additional Higgs bosons in our model. 
The pair annihilation process of ip consists of the annihilation into SM particles via Higgs- 
mediated s-channel processes and into Higgs bosons via s, t, and u-channels. The dominant 
final states of the SM particles are bb, ti, W + W~ , and ZZ. The total cross section of the 
annihilation process is given by 



av re i 



(gs sin cos I 
167r 



4m\ ' 



(3.1) 



m 



2 \2 
hi) 



+ 



1 



m 



h 2 



+ mf r 



hi hi 



2(5 - m 2 hi )(s -m\ 2 ) + 2m hl m h2 T hl T h2 



((5 — m 
2 



2 \2 



vo J 



+ 2 



m 



hi 

2s ^ 1 - 

2 



w 



2 + 



4ml 
5 



3/2 



•3 + 



■2s 1- 



4mi 



3/2 



2m 



2 ^2 
w> 



4m^ 



1 



4m 2 w 



+ 2 



m. 



vo 



2 + 



(s — 2m 2 z 
4m z 



2 \2 



4m| 



+ ^2 crh.h, + ^2 Vhihjhk, 

i,j=l,2 i,j,k=l,2 

where is the decay width of hi for i = 1, 2, and cr^hj, &hihjh k are the annihilation 
sections of tpip into hihj or hihjh^ with i,j, k = 1, 2. Here, a/s denotes the center of 
energy. The thermal average of the cross section over s is given by 



(VannVrel) ~ g . TK 2 { /T) J 



[•CO 

I ds a ann .(s) (5 - 4ml) 

J Am?, 



T 



(3.2) 
(3.3) 
(3.4) 

(3.5) 

(3.6) 
(3.7) 

cross 
mass 

(3.8) 



where K\ 2 are the modified Bessel functions. The evolution of the number density of the 
singlet fermion is described by the Boltzmann equation in terms of {(7 ann v re i) and the 
equilibrium number density of ip 



dn^ 
~dT 



+ 3Hn 



Vrel) 



(3.9) 



where H is the Hubble parameter and is the equilibrium number density of tp. After 
the freeze out of the annihilation processes, the actual number of ip per comoving volume 
becomes constant and the present relic density = m^n^ is determined. The freeze-out 
condition gives the thermal relic density in terms of the thermal average of the annihilation 
cross section. 



fi,/,/j 



(1.07 x 10 9 )z F 



fg^M p i(GeV){a ann v re i) 



(3.10) 



where g* counts the effective degrees of freedom of the relativistic quantities in equilibrium. 
The inverse freeze-out temperature xf = m^/Tp is determined by the iterative equation 



XF 



2g*x F 



\G ann.Vrel) 



(3.11) 




20 30 40 50 60 70 80 90 100 

"V (GeV) 

Figure 1: Allowed parameter set of {m^,gs) with mh 1 — 90 GeV (±1%) and m/j 2 = 500 GeV 
(±12%). The allowed region by LEP2 data is denoted as orange region. 

We investigate the allowed model parameter space, which provide thermal relic density 
consistent with the WMAP observation. In addition to and gs, we have six more 
undetermined parameters in the scalar potential: Aq, Ai, A2, A3, A4, xq, which determine 
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Figure 3: Allowed parameter set of (m^,gs) with = 120 GeV (±4%) and m^ 2 = 100 GeV 
(±1%). The allowed region by LEP2 data is denoted as orange region. 



the Higgs boson masses (m^, mh 2 ), mixing angle (9), and triple and quartic self couplings 
of Higgs bosons. Here we have to study a large volume of the multidimensional parameter 
space. For clarity of the presentation of our result, we fix the Higgs masses, and m/j 2 
within some ranges, while allowing the other parameters such as Higgs mixing angle and 
self couplings vary freely. Our parameter sets should satisfy several physical conditions. 
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We demand that i) the potential is bounded from below, ii) the electroweak symmetry 
breaking is viable, and Hi) all couplings keep the perturbativity. 

Fig.|l] shows the allowed parameter set by the measured relic abundance with m/ ll = 90 
GeV (±1%) and = 500 GeV (±12%). The valley in Fig.|l] implies the resonant region 
of hi exchange for DM pair annihilation, where 2m^ ~ mhi- The coupling constant gs 
should be small in that region in order to compensate the enhancement of the cross section 
due to the Higgs resonance effect. A step appears when ~ 80GeV, which denotes that 
the annihilation channel 'ipip — > W~W + / Z Z opens as exceeds the W and Z boson 
masses. 

The current experimental bound on Higgs mass can have a significant impact on the 
allowed parameter space. Note that the LEP2 bound of the SM Higgs boson mass gets 
weaker in our model since the SM-like Higgs couplings are modified, and therefore = 90 
GeV might be allowed depending on the parameter set. The promising channel to produce 
a neutral Higgs boson at LEP is the Higgs-strahlung process, e~e + — > Zh. A lower bound 
on the mass of the SM Higgs boson has been established to be 114.4 GeV at 95 % confidence 
level [J20y . In our model, the Higgs mixing alters the hiZZ couplings and therefore the cross 
sections of the Higgs-strahlung processes. Furthermore, Higgses can invisibly decay to a 
pair of singlet fermions. Thus the SM bound should be modified accordingly. We consider 
the parameters defined by 



/ \ 2 -pSM 

p 2 I Oh.ZZ \ L _hi 

Si I <?M I t-i.<?M , T-i/ 



(3.12) 



hi 

where are the widths of hi decays into the SM particles. Assuming the non-standard 
models, the lower bound on the Higgs mass is represented by the upper bound of £ 2 , which 
is shown in Ref. pQ] . In our analysis, we impose £ 2 < 0.1 as a conservative bound for 
m\ li = 90 GeV (and £ 2 < 0.3 for = 100 GeV for later use). Since a new fermion exists 
in our model, the definition of £ 2 includes the decay width for invisible decay channel, 
hi — > ipip. On the figure, the orange points denote the region which satisfy the Higgs 
mass bound. The allowed region appears when itu < m/ ll /2, due to large contribution 
from the invisible Higgs decay. If > m^/2, the decay channel h% — > vpip is closed and 
£i = (cos#) 2 is always larger than 0.1 so that the corresponding region is excluded by the 
experimental results of LEP2. 

Allowed parameter set with = 120 GeV and = 500 GeV are shown in Fig.|2[ 
This parameter set always satisfy the LEP2 constraints on the Higgs boson mass, because 
the Higgs masses are lager than the current experimental bound. 

Fig.||| shows the allowed parameter set for m/ ll = 120 GeV and m^ 2 = 100 GeV. Here 
m^j and mh 2 are comparable and there are two resonant regions corresponding to h\ and 
h% resonances. One can notice that the most of parameter points also satisfy the Higgs 
mass bound. 

We now briefly comment on possible LHC phenomenology of our model. One obvious 
difference of the model from SM is that we have two neutral Higgs bosons which have 
smaller couplings to SM particles, compared to SM Higgs case. Other characteristic of 
the model is that the Higgs bosons can decay invisibly to a pair of singlet fermions, if 
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Figure 4: The ratio of T^ 1 to the total decay width of hi with m hl = 120 GeV (±1%) and 
m h2 = 500 GeV (±12%). 

kinematically allowed. Fig. ^| shows the ratio of to the total decay width of hi, for 
= 120 GeV and rrih 2 = 500 GeV. We have very large invisible branching ratios when 
the mass of singlet fermion is less than half of ■ Such a large invisible Higgs decay may 
be observed at the CERN LHC pi pj. 



4. Direct detection 

There are several experiments to detect the WIMP directly through the elastic scattering 
of the WIMP on the target nuclei [p3| , |24| . The effective Lagrangian describing the elastic 
scattering of the WIMP and a nucleon is given by 



where the coupling constant f p is given by 



fp,r 



111 



EAp,n)Uq_ _2_ ( p>n ) \ - 



a a 



q=u,d,s 



q=c,b,t 



(4.1) 



(4.2) 



with the matrix elements rn^ pn ^f!f^ = (p,n\m g qq\p,n) for q = u,d,s and f^g^ = 1 



YJ 9=n d s f^ n \ The numerical values of the hadronic matrix elements f^ n> are determined 
in Ref. 



c (p,n) 



&l = 0.020 ± 0.004, f$l = 0.026 ± 0.005, /J?' = 0.118 ± 0.062 



r(p) 



(p) 



(4.3) 



and 



= 0.014 ± 0.003, = 0.036 ± 0.008, /$J = 0.118 ± 0.062. (4.4) 
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Figure 5: Predictions of the elastic scattering cross section <j(ipp — ► tpp) with respect to with 
m/ tl = 90 GeV (±1%) and m ha = 500 GeV (±12%). The red line indicates the CDMS bound, the 
purple line the Xenon bound, and the green line the up-coming super CDMS bound. The allowed 
region by LEP2 data is denoted as orange region. 



Since ftf^ dominantly contributes to /(p >n ), we let f p « f n . The effective coupling 
constant a q is defined by the spin-independent four fermion interaction of the quarks and 
the dark matter fermion in this model. The effective Lagrangian is given by 



-■int 



(4.5) 



where a q is derived by the Higgs exchange t-channel diagram to be determined by 



gs sin v cos um g 
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The elastic scattering cross section is obtained from the effective Lagrangian ( [4.1| ) 

4Af? r „ . , . AM? A 2 



(4.6) 
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TT 
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(4.7) 



where M r is the reduced mass defined by 1/M T = l/m^ + l/m nuc i e i. For the convenience of 
being compared with the experiments, we obtain the cross section with the single nucleon 
given by 



a(ipp — ► ipp) 



4wg 2 

TT /p ' 



(4.8) 



where m r is the reduced mass given by l/m r = + l/m p . 
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Figure 6: Predictions of the elastic scattering cross section a(ipp — ► tpp) with respect to with 
m hl = 120 GeV (±1%) and m h2 = 500 GeV (±12%). The red line indicates the CDMS bound, the 
purple line the Xenon bound, and the green line the up-coming super CDMS bound. The allowed 
region by LEP2 data is denoted as orange region. 

The prediction of the elastic scattering cross sections with the allowed parameter set as 
in Figjl| is depicted in Fig|| Allowed parameter sets by the LEP2 data are again denoted 
by orange points on the figure. The allowed region is entirely excluded by the current 
experiments including CDMS [^J and Xenon [28]. 

The cross sections with the allowed parameter set of Fig.|| and Fig.|| are also depicted in 
Fig.|| and Fig.|7[ respectively. The resonance region represented by a valley is not excluded 
by the experiments up to date. The large region in Fig.|| is also still allowed. We 
expect that the super CDMS experiment in the future will probe the most region of the 
parameter set for the singlet fermionic dark matter. 
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5. Conclusion 

We propose a renormalizable model with a fermionic cold dark matter. A minimal hidden 
sector consisting of a SM gauge singlet Dirac fermion and a real singlet scalar is introduced. 
We show that the singlet fermion can be a candidate of the cold dark matter which explain 
the relic abundance measured by WMAP. The constraints on the masses and couplings 
at LEP2 are included and the elastic scattering cross sections for the direct detection are 
predicted. We find that most region of the parameter set will be probed by the direct 
detection through elastic scatterings of the DM with nuclei in the near future. 
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